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Molecular dynamics simulations are used to study the occupancy and flow of water through nanotubes com¬ 
prised of hydrophobic and hydrophilic atoms, which are arranged on a honeycomb lattice to mimic function¬ 
alized carbon nanotubes (CNTs). We consider single-file motion of TIP3P water through narrow channels 
of (6,6) CNTs with varying fractions (/) of hydrophilic atoms. Various arrangements of hydrophilic atoms 
are used to create heterogeneous nanotubes with separate hydrophobic/hydrophilic domains along the tube 
as well as random mixtures of the two types of atoms. The water occupancy inside the nanotube channel 
is found to vary non-linearly as a function of /, and a small fraction of hydrophilic atoms (/ « 0.4) are 
sufficient to induce spontaneous and continuous filling of the nanotube. Interestingly, the average number 
of water molecules inside the channel and water flux through the nanotube are less sensitive to the specific 
arrangement of hydrophilic atoms than to the fraction, /. Two different regimes are observed for the water 
flux dependence on / - an approximately linear increase in flux as a function of / for / < 0.4, and almost no 
change in flux for higher / values, similar to the change in water occupancy. We are able to define an effective 
interaction strength between nanotube atoms and water’s oxygen, based on a linear combination of interaction 
strengths between hydrophobic and hydrophilic nanotube atoms and water, that can quantitatively capture 
the observed behavior. 
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I. INTRODUCTION 

Fluid flow through nanoscopic pores of carbon nan¬ 
otubes has received significant scientific attention ow¬ 
ing to the potential applications that may result from 
the use of nanotubes in various technologies 1 ®. Exam¬ 
ples include nanotube-based membranes for water de¬ 
salination and purification 5 nanoscale transport for 
lab-on-a-chip device^ and integrated circuit^. At the 
fundamental level, fluid flow through nanotubes is en¬ 
dowed with various interestin g ph enomena - breakdown 
of continuum hydrodynamics 9 23, temperature induced 
hydrophobic-hydrophilic transition^®, novel phase be¬ 
havior 1 -^, etc. Moreover, nanotubes can be used as sur¬ 
rogates to study other complex flow systems such as bio¬ 
logical channels like aquaporinS 13 . Therein, carbon nan¬ 
otubes can be used as a simple prototype devoid of com¬ 
plexity present in the biological channels, but nonetheless 
capturing some important physical detailS 14 ! 

Early molecular dynamics simulations by Hummer et 
al. showed that water can spontaneously fill the cavity of 
a (6,6) carbon nanotube and stay filled during the entire 
simulation time of 66 ns 15 . They observed a pulse-like in¬ 
termittent transmission of water molecules through the 
nanotube with similar rates as biological channels. The 
transition between empty and filled nanotubes, realized 
by tuning interactions between nanotube atoms and wa¬ 
ter®, can be described as a two-state procesd®. Several 
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theoretical models have been developed that can describe 
water permeation through the nanotube channel at equi¬ 
librium^ as well as in the presence of a driving force such 
as pressure or chemical potential difference 19 . 

While simulations on pristine CNTs have provided 
invaluable fundamental insight into water flow therein, 
functionalization of carbon nanotubes 20 is often either 
a prerequisite for their processing, e.g., the use of sur¬ 
factant for solution dispersion®, or desirable to achieve 
certain functionality such as sensitivity and selectivity in 
molecular detection®. Moreover, nanotubes with similar 
structural features as CNTs, but different atomic com¬ 
position such as Boron-Nitride 23 ^, carbon doped Boron- 
Nit rid^S, and various inorganic material^®, can now be 
synthesized routinely in the laboratory. In the context 
of CNT-surrogate study of flow through biological chan¬ 
nels, nanotube-water interactions can vary depending on 
the hydrophobicity of the amino acids involved. In all of 
these cases, as the interactions between nanotube atoms 
and water will change depending on the nanotube func¬ 
tionalization and composition, the water filling and per¬ 
meation behavior will also be different. While the under¬ 
standing of water flow through functionalized nanotubes 
so far is quite limited, simple models that incorporate 
basic details of nanotube functionalization may enable 
important progress in understanding fluid flow through 
nanotubes of a broader range of composition and func¬ 
tion. 

In this paper, we use molecular dynamics simulations 
to study water flow through model tubular channels in 
which atoms are arranged on a honeycomb lattice, and 
two different atom types are used to populate the lattice 
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sites, resulting in a binary heterogeneous nanotube. The 
two types of nanotube atoms differ only in their Lennard- 
Jones parameters with respect to water, and correspond 
to sp 2 -hybridized carbon C (labeled hydrophilic due to 
continued filling of the nanotube with water) and re¬ 
duced carbon-water interaction strength C' (termed hy¬ 
drophobic as the nanotube is mostly empty but under¬ 
goes transition between empty and filled states) as used 
by Hummer et al 15 . We note that our classification of 
nonpolar sp 2 -hybridized atoms as hydrophilic is differ¬ 
ent from previous work where the nanotube made up of 
such atoms was actually characterized as hydrophobic 15 . 
We have based our classification on the water droplet 
contact angle on a flat unrolled nanotube surface; for C 
atoms the contact angle is less than 90° (hydrophilic) and 
for C' atoms the contact angle if slightly above 90° (hy¬ 
drophobic) 16 . We first vary the fraction, /, of hydrophilic 
nanotube atoms, separated axially from the hydrophobic 
nanotube atoms, to study its impact on water filling as 
well as water flow through the nanotube. In addition, for 
a given /, we also study the influence of nanoscale pat¬ 
terning on water flow by arranging hydrophilic nanotube 
atoms in various patterns. Finally, we propose a simple 
mean-field model, in which the heterogeneous nanotube 
is represented as a homogeneous nanotube with uniform 
interactions between the nanotube atoms and water. We 
show that the predictions of this simple model are in very 
good agreement with the detailed simulation data. 


II. MODEL AND SIMULATION DETAILS 


To study the flow of water through functionalized nan¬ 
otubes, a single (6,6) nanotube, with diameter of 0.81 nm 
and length 1.34 nm, is solvated in a periodic TIP3P water 
box with 1025 molecules. The nanotube is free to trans¬ 
late and rotate througout the simulation time. The simu¬ 
lations are performed in an NPT ensemble using molecu¬ 
lar simulation package GROMACS (version 4.0.7) for 100 
ns with a timestep of 2 fs. The temperature is maintained 
at 300 K using the Nose-Hoover thermostat, and the pres¬ 
sure is maintained at 1 bar using the Parrinello-Rahman 
barostat with a time constant of 5 ps. We use the particle 
mesh Ewald (PME) method 26 for long-range electrostatic 
interactions with a short-range cutoff value of 0.9 nm. 
For short-range Lennard-Jones (LJ) interactions, we use 
a cutoff value of 1.4 nm. The well-depths of the Lennard- 
Jones potential between the nanotube atoms and water’s 
oxygen atom are 0.4784 kJ/mol (hydrophillic) and 0.2703 
kJ/mol (hydrophobic), with the corresponding LJ diam¬ 
eters being 0.3275 nm (hydrophillic) and 0.3414 nm (hy¬ 
drophobic). The angle and dihedral potentials are also 
applied to the nanotube atoms as done previously 15 . 


III. RESULTS AND DISCUSSION 

To understand the dependence of water flow through 
cylindrical channels with both hydrophobic and hy¬ 
drophilic building blocks, we use carbon nanotube like 
short cylindrical channels made up of hydrophobic and 
hydrophilic atoms. The fraction of hydrophilic atoms (/) 
is varied to tune the relative proportions of these atoms 
in order to change the water occupancy and flow char¬ 
acteristics of the resulting nanotubes. Previously, Hum¬ 
mer et al have shown that fully hydrophilic nanotubes 
(/ = 1) of the same size stay filled with an average occu¬ 
pancy of about 5 water molecules, and conduct about 17 
water molecules per nanosecond 15 . On the other hand, 
fully hydrophobic nanotubes (/ = 0) can transition be¬ 
tween filled and empty states with an apparent two-state 
kinetics 17 . Melillo et al performed a systematic varia¬ 
tion of interaction strength between nanotube atoms and 
water’s oxygen to show that water occupancy and flow 
is sensitive to chan ges in the interaction strength within 
only a narrow ranged For low attraction strengths, the 
nanotubes stay empty with no flow, whereas for high 
attraction strengths, the nanotubes are filled and con¬ 
duct water with similar rates. For intermediate attrac¬ 
tion strengths, a sharp increase in water occupancy and 
flow was observed, which was related to the change in 
wetting characteristics based on simulations carried out 
on the unrolled nanotube surface. 

Recent work from Debenedetti and co-workers has 
highlighted the importance of surface heterogeneity and 
patterning on nanoscale hydrophobicity as it relates to 
the behavior of confined water^^. Some of the com¬ 
plexities associated with water flow through heteroge¬ 
neous channels have been addressed using charged nan¬ 
otube atormf®, asymmetrically positioned charged^, and 
Y-shaped nanochannels 30 . For nanotube surfaces com¬ 
prised of hydrophilic and hydrophobic atoms, the sen¬ 
sitivity of water flow to the corresponding surface het¬ 
erogeneity, the fraction of hydrophilic atoms, and the 
actual patterning of these atoms into hydrophobic and 
hydrophilic domains, is much less clear. 

Figure [l] (left panel) shows the number of water 
molecules occupying the nanotube as a function of time 
for varying fractions of hydrophilic atoms, / = 0, 0.25, 
0.5, 0.75, and 1 (from top to bottom). In this case, hy¬ 
drophilic and hydrophobic atoms are distributed along 
the nanotube axis termed ‘axial arrangement’ to dis¬ 
tinguish it from other arrangements considered in this 
work. As expected on the basis of earlier findings, the 
fully hydrophobic nanotube (/ = 0) undergoes transi¬ 
tions between empty and filled states, whereas the fully 
hydrophilic nanotube (/ = 1) stays filled during the en¬ 
tire simulation time. Interestingly, we find that just a 
small fraction of hydrophilic atoms (/ = 0.25), localized 
at the nanotube tip, is sufficient to significantly increase 
the number of water molecules occupying the nanotube. 
For / > 0.4, the nanotube stays completely filled. 

To further quantify water occupancy as a function of 
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FIG. 1. Water occupancy. (Left panel) Number of water molecules N for varying fraction of hydrophilic nanotube atoms / are 
shown as a function of time. (Middle panel) Nanotube snapshots are shown, with hydrophobic and hydrophilic atoms labeled 
in red and blue, respectively. (Right panel) Probability of observing N water molecules inside the nanotube P(N) is shown for 
various values of /. The inset shows average number of water molecules as a function of / for three possible arrangements of 
hydrophilic atoms. 


/, we calculate the distribution P(7V), which is the prob¬ 
ability of finding N water molecules inside the nanotube. 
The P(N) for several representative values of / are shown 
in Fig. [I] (right panel). The probability of observing an 
empty nanotube, P(0), decreases dramatically with in¬ 
creasing number of hydrophilic nanotube atoms (i.e., in¬ 
creasing / to 0.5). This precipitous decrease in P(0) is 
complemented by a similarly dramatic rise in the proba¬ 
bility of observing N = 5 or 6 water molecules, whereas 
relatively little change in P(N) is observed with subse¬ 
quent increase in / from / = 0.75to/ = l. As shown 
in the inset of Fig. [l] (right panel), the average number 
of water molecules inside the nanotube < N > changes 
approximately linearly as a function of / for / < 0.4 and 
then plateaus at a value found for the fully hydrophilic 
case. The so-called ’axial arrangement’ is the same as 
depicted in Fig. [TJ which is obtained by separating con¬ 
tinuous domains of hydrophilic and hydrophobic atoms 
along the nanotube axis. The so-called ‘random arrange¬ 
ment’ is obtained by placing the hydrophilic atoms ran¬ 
domly on the nanotube lattice (lattice sites numbered 
from 1 to 144), with two different cases considered here 
by changing the seed of the random number generator. 

In addition to water occupancy, water transport and 
its potential sensitivity to tube properties, are important 
to quantify. To probe the effects of changing the fraction 
of hydrophilic atoms on water flux through the nanotube, 
we calculate the number of water molecules that perme¬ 
ate the nanotube by entering either end, but leaving from 
the opposite end over the total simulation time. Water 
flow is defined as the number of water molecules trans¬ 
mitted per nanosecond through the nanotube. Figure [2] 
shows water flow as a function of / for three possible 
arrangements of hydrophilic nanotube atoms. For refer¬ 
ence, we show the value of water flow obtained by Hum¬ 
mer et al 15 for the fully hydrophilic nanotube (/ = 1). 

First, focusing on the data for the axial arrangement 
of hydrophilic and hydrophobic atoms in Fig. |2j we ob- 
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FIG. 2. Water flow rate. Number of water molecules perme¬ 
ating through the nanotube, which enter either side and leave 
from the opposite side, per unit time, are shown as a function 
of fraction of hydrophilic atoms /. The symbols connected by 
dashed lines are the simulation data and the solid black line 
is the linear fit for / < 0.4 for the ‘axial arrangement’ data. 

serve an approximately linear dependence of water flow 
on / for / < 0.4. Outside this regime, the water flow is 
relatively insensitive to the increase in /, and a slight 
decrease in water flux is observed with increasing /. 
The linear dependence of water flow as a function of 
/ for / < 0.4 likely appears due to the concomitant 
linear relationship between water flow and occupancy, 
shown for the same condition in Fig. |T] (inset). A sim¬ 
ilar change in water occupancy and flow was reported 
previously for a spatially uniform change in nanotube- 
water attraction strength over low to moderate interac¬ 
tion strengths 1 ® There, an increase in the nanotube- 
water attraction strength was shown to monotonically 
increase water’s residence time inside the nanotube. For 
higher attraction strengths, this increase in the residence 
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FIG. 3. Number of water molecules permeating through the 
nanotube, which enter either side and leave from the opposite 
side, per unit time are shown for various arrangements of 
hydrophilic (blue) and hydrophobic atoms (red) for / = 0.5 
as depicted in the nanotube snapshots. 


time led to a decrease in water flow, thereby deeming 
both the average water occupancy and the residence time 
important for capturing the change in water flow. In the 
case of heterogeneous nanotubes, an increase in the frac¬ 
tion of hydrophilic atoms above / = 0.4 leads to a slight 
increase in water occupancy (Fig. [I), but insensitive or 
subtly decreasing water flow (Fig. 2]). Therefore, the 
water flow characteristics appear to be influenced sim¬ 
ilarly by a heterogeneous distribution of hydrophilic and 
hydrophobic atoms as compared to a spatially uniform 
change in nanotube-water interaction strength. 

To see if the specific axial arrangement of hydrophilic 
atoms is responsible for the observed behavior, Fig. [2] also 
shows water flow as a function of / for random arrange¬ 
ments of hydrophilic atoms. The qualitative behavior 
is quite similar in all cases, with approximately linear 
dependence for low / values and a plataeu (or slight de¬ 
crease) for higher / values. 

Though statistically less significant, a slight increase 
in water flow for a given fraction of hydrophilic atoms 
is observed in the case of the random arrangement as 
opposed to the axial arrangement. The water occupancy 
data (Fig. 1, right panel inset) also shows a similar trend, 
with the random arrangement leading to slightly higher 
water occupancy. One can explain this based on previous 
work by Giovambattista et al. where it was shown that 
hydrophilic borders surrounding nanoscale hydrophobic 
patches can significantly alter water hydration by reduc¬ 
ing the water repulsion 271 . In the ‘random arrangement’ 
case, the hydrophilic atoms are mixed with hydrophobic 
atoms, thereby enhancing water density near such sur¬ 
faces, which can lead to slightly enhanced water flow. 

Between the cases of axial and random arrangement 
of hydrophilic and hydrophobic atoms exists a range of 
cases wherein discrete domains could conceivably pattern 
the nanotube surface. To further probe how such specific 
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FIG. 4. Flow through homogeneous versus heterogeneous 
nanotubes. (Top panel) Number of water molecules permeat¬ 
ing through the nanotube, which enter either side and leave 
from the opposite side, per unit time are shown as a function 
of interaction strength between nanotube atoms and water’s 
oxygen. For nanotubes comprised of hydrophobic and hy¬ 
drophilic atoms, we use eq. 1 to estimate an effective interac¬ 
tion strength e e ff, whereas separate simulations are conducted 
for nanotubes comprised of a single type of atoms interacting 
with water’s oxygen with cnt-ow- (Bottom panel) Aver¬ 
age number of water molecules as a function of interaction 
strength between nanotube atoms and water’s oxygen. 


arrangement of hydrophilic atoms (e.g., on the ends of the 
nanotube) may modulate water flow, we simulate several 
possible arrangements as shown in Fig. [3j all for / = 0.5. 
Within statistical uncertainty of the measured water flow 
rate, the different heterogeneities studied produce similar 
results. As the water flow through these short nanotubes 
is expected to be dominated by the barriers at the en¬ 
try and exit of the nanotubc 5 i the observed insensitivity 
of water flow to the actual arrangement of hydrophilic 
atoms is a bit surprising. While this effect may be tuned 
by increasing nanotube length and for low values of /, 
this requires future work employing either a low-cost but 
accurate single-site water model 31 or dipole lattice model 
of water 32 l 

The qualitative similarity we observe here in water flow 
through heterogenous nanotubes made up of two types 
of atoms with previous work on homogeneous scaling of 
nanotube-water interaction strength, begs the question 
about whether these two cases may be related in a more 
quantitative way. Due to the observed insensitivity of the 
water occupancy and flow on the actual arrangement of 
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the hydrophilic atoms, we hypothesize that one can define 
an effective interaction strength e e ff between nanotube 
atoms and water’s oxygen as, 

^eff ^philic f T ^phobic(1 /)? (1) 

where e p hm c and e p hobic correspond to the LJ inter¬ 
action parameters between nanotube atoms and water’s 
oxygen for hydrophilic and hydrophobic atoms, respec¬ 
tively. Such a linear interpolation is also justified based 
on the small differences in the interaction potential curves 
(LJ potential) between the two cased^. To test if such 
a simple scaling relationship works over a broad range, 
we calculate water flow through a nanotube and average 
water occupancy for various nanotube-water interaction 
strengths cnt-ow* The comparison of this data with wa¬ 
ter flow data from Fig. [2] and water occupancy data from 
Fig. [l] (right panel inset) is shown in Fig. [4j Overall, the 
change in water flow and < N > as a function of uni¬ 
form change in nanotube-water interaction strength, as 
predicted by detailed simulations, is in remarkably good 
agreement with the prediction based on the mean-field 
type description in Eq. 1 for heterogeneous nanotubes. 
This quantitative agreement in Fig. [4] suggests that the 
fluid flow through functionalized nanotubes may be a 
simple function of a parameter describing the affinity be¬ 
tween the fluid molecules and the nanotube surface. The 
apparent insensitivity to the distribution of the two atom 
types may be a manifestation of the single-file nature of 
water flow through a (6,6) nanotube. This highlights 
the question for future study regarding whether devia¬ 
tions from a mean-field description become significant 
for larger diameter nanotubes or not. In addition, it will 
be useful for future studies to elucidate any sensitivity of 
the observed behavior to nanotube length. 

IV. CONCLUSIONS 

We have studied the water filling and transport 
through a simple model of functionalized carbon nan¬ 
otubes. The fraction of hydrophilic atoms, /, are varied 
to change the nanotube from fully hydrophobic (/ = 0) 
to fully hydrophilic (/ = 1). Two different regimes are 
identified for water occupancy and flow as a function of 
/. For low / values (/ < 0.4), there is an approximately 
linear dependence on / of the average number of water 
molecules inside the nanotube as well as the number of 
water molecules transmitting through the nanotube. For 
higher / values, a slight increase in water occupancy is 
observed, whereas a slight decrease in water flow is pre¬ 
dicted. For most of the systems considered in this work, 
the observed water flow properties are relatively insen¬ 
sitive to the precise arrangement of hydrophilic atoms. 
A simple mean-field type expression for effective inter¬ 
action strength between nanotube atoms and water can 
be used to quantitatively capture flow characteristics of 
water through heterogeneous nanotubes. These findings 


suggest that control over the character (i.e., interactions 
with water) of the nanotube functionality may be more 
important than the specific spatial arrangement of that 
function for controlling water flow therein. This simple 
design principle should enable prediction of flow char¬ 
acteristics in synthetic and functionalized nanotubes be¬ 
yond pristine CNT and help guide rational design of tar¬ 
get functionality. 
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